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Abstract 

In this paper a search for associated charm production both in neutral and charged current ^-nucleus 
interactions is presented. The improvement of automatic scanning systems in the CHORUS experi- 
ment allows an efficient search to be performed in emulsion for short-lived particles. Hence a search 
for rare processes, like the associated charm production, becomes possible through the observa- 
tion of the double charm-decay topology with a very low background. About 130,000 v interac- 
tions located in the emulsion target have been analysed. Three events with two charm decays have 
been observed in the neutral-current sample with an estimated background of 0.18±0.05. The rel- 
ative rate of the associated charm cross-section in deep inelastic v interactions, a(ccv)la^^ — 
(3.62^2'42( s t a t) ± 0.54(syst)) x 10~ 3 has been measured. One event with two charm decays has 
been observed in charged-current interactions with an estimated background of 0.18±0.06 and 
the upper limit on associated charm production in charged-current interactions at 90% C.L. has been 
found to be a(cc[j,~)/ ace < 9.69 x 10~ 4 . 
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1 Introduction 

The production of two charmed particles in neutrino-nucleon scattering is a very rare process and 
therefore difficult to measure. The charm pair originates from two different processes: one, the so-called 
boson-gluon fusion mechanism, is possible only in neutral-current interactions, while the other, gluon 
bremsstrahlung, occurs both in neutral and charged current (NC and CC) interactions. The Feynman 
diagrams of both process are shown in Figure Q] 

In early neutrino experiments at the CERN SPS and Fermilab, measurements of production rates 
of prompt trimuons [T] and like-sign dimuons (0,0 and 0) were found to be higher than expected 0. 
In the past the origin of these events was assumed to be the associated charm production in CC neutrino 
interactions with subsequent muonic decay of the charmed particle(s): 

v^N -» iTccX (1) 

where the charm-anticharm pair is produced by a gluon emitted via bremsstrahlung from a light partem. 
More recently, an analysis by CHORUS [6] suggested a significant contribution from decays of 

light vector mesons. In a previous search, based on the visual observation of the charmed hadron decays, 
the CHORUS collaboration has published the first evidence for this process [7] through the observation 
of one event. 

For the NC process, the E531 experiment reported the observation of one event with two neu- 
tral decay topology. Using this event and correcting for the efficiency, a relative rate was quoted as 

l^f^^J = ( L3 ±li) x 10_3 [ § 1 at the average neutrino beam energy of 22 GeV. No background 
estimation was given. 

A measurement of this rate was recently made by the NuTeV experiment [9 ] through the detection 
of events with wrong sign muon final states. The analysis technique consists in comparing the visible 
inelasticity, y^- ls = Euad/ (Ehad + E^), measured in the and data samples where the leading 
muon had the wrong sign of charge, to a Monte Carlo simulation containing all known conventional 
wrong sign muon sources and a possible NC charm signal. The NC charm signal peaks at large values 
of y v - ls because the muon from the decay of a heavy flavor hadron is usually much less energetic than 
the hadron shower produced in the NC interaction. After consideration of all background sources, a 
production cross section ^^r^TT = (6.4*!' g) x 10~ 3 at the average neutrino beam energy of 154 
GeV was obtained. More details both on the theoretical and experimental aspects of associated charm 
production may be found in Ref . iflOl [TTIl and references in them. 

All the results, available to date, with the exception those from E531, were extracted from multi- 
muon samples with a large background of non prompt muons. In the CHORUS experiment a different 
search for associated charm production in NC and CC interactions was carried out in nuclear emul- 
sion, is practically background free and is presented in this paper. The search is based on the visual 
observation of the primary vertex and the two decay vertices, taking advantage of the sub-micrometric 
position resolution of the nuclear emulsion. We report the observation of four events with two charm de- 
cays: three in NC and one in CC interactions. Based on the observed events, taking care of the estimated 
background and correcting for the detection efficiencies, we obtain a production rate of associated charm 
in NC interactions and an upper limit for the CC process. 

2 The experimental apparatus 

The CHORUS detector is a hybrid setup that combines a nuclear emulsion target with various 
electronic detectors, as described in Ref. [12J. 

A three-dimensional reconstruction of the decays of short-lived particles such as the r lepton and 
charmed hadrons can be performed in nuclear emulsions. The emulsion target is segmented into four 
stacks with an overall mass of 770 kg. Each of the stacks consists of eight modules of 36 plates of 36 cm 
x 72 cm size. Each plate consists of a 90 fim plastic support coated on both sides with a 350 fim emulsion 
layer 1 13 ]. Each stack is followed by three emulsion sheets, having a 90 fim emulsion layer on both sides 
of an 800 ^m thick plastic base, acting as an interface to a set of scintillating fibre tracker planes. The 
accuracy of the fibre tracker predictions is about 150 fim in track position and 2 mrad in angle. The 
interface sheets and the fibre trackers provide accurate particle trajectory predictions in the emulsion 




Figure 1: Feynman diagrams for boson-gluon fusion(left) and gluon bremsstrahlung processes(right). 

stack in order to locate the interaction vertex. The information of the electronic detectors has been used 
to define two data sets, the Ifj, and Ofi samples. The events belonging to the 1/i sample contain one 
reconstructed muon of negative charge. The muon identification and reconstruction is based on the muon 
spectrometer response. The 0/i sample contains events where no muon is found. It mainly consists of NC 
U/j, interactions with a contamination of mis-identified CC interactions and interactions generated by 
neutrinos other than v^. 

The emulsion scanning is performed by fully automated microscopes equipped with CCD cameras 
and a read out system called "Ultra Track Selector" [14]. In order to recognize the track segments in 
the emulsion, a series of tomographic images are taken by moving the focus at different depths in the 
emulsion thickness. The digitized images are shifted according to the predicted track angle and then 
added. The presence of aligned grains forming a track is detected as a local peak above the gray level of 
the summed image. The track finding efficiency of the track selector is higher than 98% for track angles 
up to 400 mrad with respect to normal incidence. 

The electronic detectors downstream of the emulsion target include a hadron spectrometer which 
measures the bending of charged particles in an air-core magnet, a calorimeter where the energy and 
direction of showers are measured and a muon spectrometer which measures the charge and momentum 
of muons. In this analysis, the momentum of particles other than muons in the candidate events were 
determined by measuring their multiple scattering in the emulsion target. 

3 Event samples and selection of two charm topologies 

The West Area Neutrino Facility (WANF) at CERN provided a wide band beam of 27 GeV average 
energy consisting mainly of u^, with a contamination of 5.1% v^, 0.8% v e and 0.2% v e . During the four 
years of CHORUS operation the emulsion target has been exposed to the beam for an integrated intensity 
which corresponds to 5.06 x 10 19 protons on target. The data from the electronic detectors have been 
analyzed and the set of events possibly originating in the emulsion stacks has been identified. 

In total about 130,000 v events have been located in emulsion with a procedure described in 
Ref. iPTSll . In order to identify a charm decay, a volume of 1.5 mm x 1.5 mm x6.3 mm around the located 
vertex position is scanned. The parameters of all track segments with angles below 400 mrad found in 
this volume are stored in a database. Typically, five thousand track segments are recorded per event. This 
procedure is called "Netscan" lfl6l . A detailed description of the procedure and of the algorithms used to 
reconstruct the event topology are given in Ref. lfT71[T8l . 

Track segments collected in the Netscan fiducial volume are analysed off-line to reconstruct the 
complete event topology: segments from different plates are combined into tracks; low-energy tracks are 
filtered out as well as track originating from outside the volume (mainly passing-through beam muons). 

Possible vertices are then defined, requiring that the minimum distance between tracks attached to 
a common vertex is compatible with zero within errors. If several vertices are reconstructed, the primary 
vertex is defined, for the 1/z sample, as the one to which the muon track is attached. For the 0/x sample, 
the primary vertex is chosen as the most upstream one with at least one track matching one of the tracks 



Table 1: Results of the visual inspection in the selected 0/x sample. 



Charm candidates 


Rejected events 


Topology 


Events 


Category 


Events 


V2 


145 


Low momentum 


72 


V4 


43 


Traversing tracks 


108 


CI 


80 


/i ± int. 


110 


C3 


75 


h° int. 


4 


C5 


4 


7- conversion 


36 


V2+V2 


1 


Overlay secondary vertices 


37 


C1+V2 


1 






C3+V4 


1 






Total 


350 


367 



reconstructed in the electronic detectors. Events for which the emulsion tracks do not converge to a single 
common vertex potentially contain a short-lived particle decay topology. They were selected for further 
inspection if one of the following selection criteria is satisfied fP71[T8l : 

- Two or more multi-tracks vertices are reconstructed. 

- At least one track of the primary vertex is not reconstructed up to the edge of the Netscan volume 

- A track originating in the same plate as the primary vertex has an impact parameter IP relative to 
it satisfying the two conditions: 

• IP > a/(5.0 2 + (2dx x a) 2 ) /jm where a = ^/(O.OOS 2 + (0.0194 x tan(9) 2 ) is a parame- 
terization of the angular error and dx is the distance of the vertex from the most upstream 
daughter track segment. 

• IP < 750 fim, to reject spurious tracks not related to the neutrino interaction. 

After selection, the samples of 26,621 Ofi and 99,245 events reduce to respectively 717 and 
2816 events. These are then visually inspected to confirm the decay topology. A secondary vertex is 
accepted as decay if the number of prongs is consistent with charge conservation and no evidence for 
hadronic interaction, e.g. Auger electron, nuclear recoil or break up, is observed. Decays into a single 
charged particle (so called kink) are accepted only if the angle between the parent and daughter tracks is 
greater than 50 mrad and the length of the parent track is greater than 25 /im. 

The observable decay topologies are classified according to their numbers of charged decay prod- 
ucts into CI, C3 or C5 for the charged and V2, V4 or V6 for neutral decaying particles. The result of the 
visual inspection is given in Table Q] and Table [2] for the Ofi and the Ifi samples respectively. In both the 
O/i and the 1/x samples, three events have been observed with two charm-decay topology. 

The rejected sample consists mainly of low momentum tracks for which the multiple scattering 
leads to a fake impact parameter or to a failure to reconstruct the track as traversing the plate. The re- 
maining events in the sample consist of hadronic interactions, gamma conversions and overlay secondary 
vertices, reconstructed using one or more background tracks or of vertices with a parent track not con- 
nected to the primary vertex. 

After topological confirmation of the decays, further kinematical selections are applied in order 
to eliminate the background coming from strange particle decays and from hadronic interactions. In 
case of a CI topology, the requirement > 0.25 GeV/c on the transverse momentum of the daughter 
particle with respect to the parent direction is applied. This cut also reduces significantly the background 
caused by "white kinks", i.e. single prong hadronic interaction without emission of any heavily ionizing 
particle, Auger electron or other evidence for nuclear recoil or break up. The selection criterion applied 
to V2 decays is >10 mrad, where the angle <j) measures the acoplanarity between the parent direction 
and the plane formed by the charged decay products. This cut eliminates mainly the two-body decays of 
neutral strange particles. The effects of these kinematical cuts are included in the determination of the 
Netscan efficiencies. 

Four events satisfied the above selection criteria. Their general features are listed in Tables [3] and 

SI 



Table 2: Results of the visual inspection in the selected 1/i sample. 



Charm candidates 


Rejected events 


Topology 


Events 


Category 


Events 


V2 


841 


Low momentum 


140 


V4 


230 


Traversing tracks 


149 


V6 


3 


/i ± int. 


258 


CI 


462 


h° int. 


69 


C3 


501 


7- conversion 


101 


C5 


23 


Overlay secondary vertices 


36 


C1+C3 


1 






C1+V2 


1 






V2+V4 


1 






Total 


2063 


753 



4 The candidate events in the NC sample 

In the NC sample three candidates for associated charm production satisfied the above selection 
criteria. A description of each event is given below. 

Event 8132-12312 with a V2+V2 decay topology is sketched in Figure|2]and microscope view of 
the event in the emulsion at relevant positions is shown in Figure|3] Only one minimum-ionizing (shower) 
particle is emitted from the primary vertex and is identified as a hadron in the electronic detector. In the 
same emulsion plate, a neutral particle decays into two charged particles (V2) 63 fim downstream of 
the primary vertex. Both daughter tracks are reconstructed in the electronic detector. The acoplanarity of 
parent and daughter tracks, <fi = 24.2 mrad ± 5.3 mrad rules out a two body decay. Moreover, the lower 
limit of the mass M m - m of the parent particle is measured to be 0.80 GeV (at 90% C.L.). In one plate 
downstream, a second V2 decay has been located, 977 /xm from the primary vertex. The acoplanarity 
angle <ft and the minimum mass M m ; n are measured to be 36.1 mrad ± 0.3 mrad and 0.95 GeV (at 90% 
C.L.) respectively. 

Table 3: General features of the associated charm candidates. Ns(Ns*) and Nh stand for number of shower 
(reconstructed in the electronic detector) and heavily ionizing tracks at the neutrino interaction vertex. 
E s h is the hadronic shower energy, p^ the muon momentum and #^(y,z) the angle between the muon and 
the beam directions, respectively in the horizontal and vertical planes. No correction was made for the 
energy of the two neutrinos from the charm decays. 





Event Id 


Ns (Ns*) 


Nh 


E sh (GeV) 


p M (GeV/c) 


fl M (y,zXrad) 




8132-12312 


1(1) 


2 


38.2 






NC 


7692-5575 


2(1) 





47.1 








7739-3952 


6(1) 


1 


45.9 






cc 


7904-4944 


4(2) 





42.5 


17.6 


(-0.164,0.080) 



Event 7692-5575 (with a C1+V2 decay topology): At the primary vertex there are two shower 
tracks, one of which decaying after 163 fim into a charged particle (C 1) with a 0. 179 rad kink angle. The 
daughter momentum was measured by multiple scattering to be larger than 1.82 GeV/c (at 90% C.L.). 
Its transverse momentum p^ with respect to its parent is larger than 0.32 GeV/c (at 90% C.L.), which 
rules out the hypothesis of a strange particle decay. In the same emulsion plate, 224 ftm from the primary 
vertex, a second decay with a V2 topology has been observed; the decay vertex occurs however in the 
plastic base. The acoplanarity angle is 14.9 mrad ± 1.5 mrad and M m ; n = 0.68 GeV (at 90% C.L.) These 
values are inconsistent with a strange-particle decay hypothesis. 

Event 7739-3952 (with a C3+V4 decay topology): The primary vertex has six shower tracks, none 
of them identified as a muon. One of the shower tracks decays after 426 fim into three charged particles 



Table 4: General features of the charged daughter particles in the associated charm candidates. M m - m 
stands for minimum mass at 90 % C.L. of the parent particle. 9(y,z) is the angle between the charged 
daughter particle and the beam directions, respectively in the horizontal and vertical planes, p is the mo- 
mentum of the charged daughter particle measured using multiple coulomb scattering. (The momentum 
measurement was not applied to Event 7739-3952 and failed for two daughters of V4 decay in Event 
7739-3952.) 





Event Id 


Topology 


M min (GeV) 


%,z)(rad) 


p (GeV/c) 


NC 


8132-12312 


V2 


0.80 


(0.049 -0.026) 
(-0.023,0.040) 


4.66 
1.39 






V2 


0.95 


(-0.165,0.045) 
(0.034,0.105) 


2.72 
3.33 




7692-5575 


V2 
CI 


0.68 


(0.026,-0.056) 
(0.082,0.032) 
(-0.136,0.291) 


4.11 
3.81 
1.82 




7739-3952 


V4 
C3 




(0.129,-0.072) 
(-0.005,-0.064) 
(0.083,0.042) 
(0.075,0.336) 
(0.011,-0.038) 
(-0.071,0.238) 
(-0.294,-0.376) 




cc 


7904-4944 


V4 




(0.007,-0.087) 
(-0.045,-0.001) 
(0.269,0.305) 
(0.101,-0.169) 


1.92 
2.32 






V2 


0.81 


(0.029,-0.053) 
(0.001,0.289) 


4.70 
0.67 



while the V4 decay vertex is 884 fim from the primary vertex. 

5 The candidate event in the CC sample 

Out of the three double charm candidates of Table 2, two are rejected on the basis of the p^ cut 
at their C 1 vertex. The remaining event 7904-4944 is shown in Figure 0] The primary vertex is located 
200 [im from the plate downstream surface. There are four primary shower tracks; one of them is identi- 
fied as a negative muon in the spectrometer. In the same plate, 58 fim downstream of the primary vertex, 
a neutral particle decays into a V2 topology. Both daughter tracks are reconstructed in the electronic 
detector. The acoplanarity angle <fi is 12.8 mrad ± 1.5 mrad, and the minimum mass of the parent par- 
ticle, M m i n = 0.81 GeV at 90% C.L. A second neutral particle decays into a V4 topology with a flight 
length of 761 ^m from the primary vertex. Two of four daughter tracks are reconstructed in the electronic 
detector. 

In Ref. Q, we have reported the observation of one event in CC interactions based on a subset 
of the data (about 50% of the present sample). This search was based on different selection criteria 
and was performed before the development of the Netscan technique: charged particles at the primary 
vertex were followed down manually to find decay vertices. This event has not been retrieved in the 
present analysis because the emission angle of the charm daughter is out of the angular acceptance of the 
Netscan procedure. 

6 Evaluation of the reconstruction efficiency 

Large samples of neutrino interactions were generated with a neutrino energy distribution ac- 
cording to the beam spectrum using the HERWIG event generator lfT9l . with the leading-order partem 
distribution functions of MRST ll20l . This generator produces associated charm-production events in 




both neutral-current and charged-current interactions. Efficiencies and backgrounds were evaluated with 
a simulation of the detector based on GEANT3 ll2"Tll . 

The simulated response of the electronic detectors is processed through the same reconstruction 
program as the one used for the data. The tracks in emulsion and the performance of the UTS are also 
simulated, in order to evaluate the efficiency of the scanning procedure. 

In order to evaluate the Netscan procedure efficiency, one needs to reproduce realistic track den- 
sities. This was achieved by merging the emulsion data of the simulated events with real Netscan data 
which do not have a reconstructed vertex but contain tracks which stop or pass through the Netscan fidu- 
cial volume, representing the real experimental background l22l . The combined data are passed through 
the same Netscan reconstruction and selection programs as used for the real data. The average detection 
efficiencies for NC and CC sample are evaluated to be (12.8 ± 0.6)% and (8.5 ± 0.6)%, respectively. 

The ratio of the reconstruction and location efficiency of the NC cc events to that of the whole 
sample of NC events is found to be 0.496 ± 0.022. This ratio is 0.506 ± 0.018 in the CC channel. 

The lower reconstruction and location efficiency of cc events is due to the larger hadronic activity 
in the electronic detector for these events. 



Table 5: Background sources and the corresponding event yields in the 1/i sample. 



Primary interaction secondary vertex event yield 





-» ii 


"cS ± X 


Y, ± — » 1 prong 


0.002 ± 0.001 


v„N- 


-» \t 


~ch ± X 


white CI 


0.08 ±0.04 


v»N- 


* M~ 


~ch ± X 


/i ± white C3 


0.06 ± 0.04 


v»N- 


-» \i 


-ch ± X 


/i ± white C5 


0.024 ± 0.009 


v»N- 


-> fj,' 


-ch°X 


h° white V2 


0.01 ± 0.01 




+ li 


-ch°X 


h° white V4 


0.006 ± 0.006 



Total 0.18 ±0.06 



7 Background evaluation 

The search is essentially topological even if there are additional measurements of the momentum 
and acoplanarity angle to further validate the candidates. In Table [5] we list all relevant background 
sources for double charm production in CC interactions. 

A double charm event is topologically indistinguishable from a single charm event with a primary 
non-charmed hadron which undergoes either a decay or an interaction without any visible recoil of the 
nuclei. Among the strange hadrons produced in a single charm event, the S 1 * 1 has the lowest cr (a few cm) 
and is relatively abundant (about 1%). Owing to longer decay lengths there are negligible contributions 
from decays of other non-charmed hadrons like kaons and pions. 




Figure 3: Microscope views of Event 8132-12312 taken in the emulsion. Since the neutrino beam is 
perpendicular to the emulsion target, the trajectory of a charged particle appears as a 'dot'. The top- 
left picture shows the neutrino interaction vertex. A charged particle which is visible several microns 
downstream of the vertex emerges from the neutrino vertex. The top-center picture shows the decay 
vertex of a neutral charmed particle decaying into two charged particles (V2) 63 fim downstream of the 
primary vertex. The trajectory of the daughter particles is visible at bottom-center picture. A second V2 
decay topology is found 977 /im downstream of the primary vertex as shown in the top-right picture. 
In the bottom-right picture the trajectory of daughter particles from the second neutral decay becomes 
visible after several microns downstream of the decay vertex. 

Decays of strange neutral particles like A and do not give a sizeable contribution owing to the 
acoplanarity cut. 

The interaction length for white interactions is very large, ranging from a few metres for single 
prong interactions to several hundred metres for multi-prong interactions. A simulation of these processes 
using FLUKA [ 23 ] has been carried out and the interaction lengths for the different processes have been 
obtained. In the course of the CHORUS analysis, 243 m of charged hadron tracks have been followed 
and 26 white interactions iPTSI have been collected, which allows a cross-check of the simulation at the 
35% level. An overall background of 0.18 ± 0.06 events is expected in CC interactions, mainly from 
white interactions. 

In Table [6] are listed the background sources and the corresponding event yield in the Ofi sample. 
An overall background of 0.18 ± 0.05 events is estimated. The main background source for the double 
charm production in NC interaction comes from charged lepton misidentification in CC double charm 
events. This background is estimated by scaling the candidate event in CC interactions with the CC 
contamination in Omu sample which is estimated to be about 14%. 

8 Results and conclusions 

In summary, in a sample of 99,245 lfi events, one event has decay topologies that meet the re- 
quirements for associated charm production. The total background in the CC sample is estimated to be 
0.18 ± 0.06 events. In the sample of 26,621 0/j events, three events showed decay topologies that are 




Figure 4: Sketch of CC event (7904-4944) with a V2+V4 topology (the shaded area is the 90 /im plastic 
base between the two 350 emulsion layers). 



Table 6: Background sources and the corresponding event yields in the Ofi sample where (jl (e ) is 
misidentified. 



Primary interaction secondary vertex event yield 



vN - 


-> M ± (e ± 




S ± — > 1 prong 


0.0003 ± 0.0001 


uN - 


-> // ± (e ± )c h±X 


/i ± white CI 


0.013 ± 0.005 


uN - 


-> M ±(e± 




white C3 


0.01 1 ± 0.005 


uN - 


- f(e± 




/i ± white C5 


0.004 ± 0.001 


uN - 


-> /x±(e± 


)c/i°X 


h° white V2 


0.001 ± 0.001 


vN - 




)ch°X 


white V4 


0.0006 ± 0.0006 


vN - 




)ccX 


cc — > any prong 


0.146 ±0.051 



Total 0.18 ±0.05 



consistent with associated charm production. The total background in the NC sample is estimated to be 
0.18 ±0.05. 

In order to estimate the associated charm production rate in interactions, an additional weight 
factor needs to be applied to lfi events with > 30 GeV/c, since a small fraction of this category 
was not located and analysed. This factor was evaluated from the measured ratio 0.305 ± 0.002 of 
events with > 30 GeV/c to those with p^ < 30 GeV/c; it was found to be 1.021. In order to evaluate 
the real number of v NC we have used the measured fraction of 0.274 ± 0.005 between deep-inelastic 
NC and CC interactions. This value is smaller than what was measured by CHARM, CDHS and 
CCFR collaborations [24 ] since a correction for the non-isoscalarity of the emulsion target is applied. 
In NC interactions identification of the neutrino flavor is not possible in the experiment. Therefore, the 
NC sample contains a small fraction of u e , v e and NC interactions. The contribution of neutrino 
flavors other than to the normalization is estimated to be 3.3%. Since in CC interactions at least one 
reconstructed muon in the spectrometer is required both NC and CC interactions have a similar energy 
threshold which is about 6 GeV. 

The relative rate of NC associated charm production is given by 

ajccv) _ N£ 8 -N hgr 1 1 
o"^cf R^Ncc n oc e net ' 

where: 

- A^ g = 3 is the number of candidate events in the 0/x sample; 

- iVbg r = 0.18 ± 0.05 is the total background in NC sample; 

- i?Nc = 0.285 ± 0.005 is the effective ratio between deep-inelastic NC and CC interactions; 

cc 

- iVcc = 101,329 is the number of CC interactions; 

- £nct = 0.128 ± 0.006 is the average detection efficiency for the NC cc sample. 



- n oc = 0.21 1 ± 0.004 is the ratio of the reconstruction and location efficiency of events with cc in 
NC events to that of all CC events. 
The value obtained for this ratio normalizing to the total neutrino flux with 27 GeV average neutrino 
energy is 

-^jf- = (3.62+|^ (stat) ±0.54(syst)) x 10~ 3 . 
°nc 

The statistical error is derived using a 68% confidence interval in the unified approach for the analysis 
of small signals in the presence of background |[27l . We have accounted for a systematic uncertainty of 
15% coming from the efficiency estimation by Monte Carlo modeling. The energy threshold of associated 
charm production is significantly higher than for NC interactions. Based on the HERWIG event generator, 
an energy threshold cut of 35 GeV is applied to NC interactions. The relative rate of NC associated 
charm production is also evaluated with the energy threshold cut and found to be (7.30l4'||(stat) ± 
1.09(syst)) x 10~ 3 . The average neutrino energy above this threshold is 73 GeV, to be compared with the 
average visible energy of our candidates of 44 GeV. The neutrino spectrum for NC interactions is peaked 
at about 25 GeV and decreases very rapidly. Therefore, we expect events to occur close the threshold. 
The result is consistent with the E531 [8] and NuTeV [9] measurements. The measured cross-section is 
consistent with our measurement on NC production of J/^ ll25lPl . Given the fact that the cross-section is 
predicted to have a strong energy dependence and taking into account that the neutrino energy spectrum 
of the beam used in this experiment peaks at values where this cross-section is predicted to be low, our 
result is compatible with the prediction of the Z°-gluon fusion model ll26l . In the framework of this 
model, one obtains the relative rate of NC associated charm production as ~ 4 x 10~ 3 . 

With the observation of one event for associated charm production in CC interactions, we obtain 
for the relative rate an upper limit at 90% C.L. ll27l of 

^Ll < 9 . 69 x io-4, 
C"CC 

normalizing to the total neutrino flux. Due to different energy thresholds in single and associated charm 
production, the production cross-section relative to CC interactions is estimated with the energy cut of 
35 GeV. An upper limit for the relative rate with the energy threshold cut is found to be 2.24 x 10~ 3 . The 
average neutrino energy above this threshold is 73 GeV. Since the topology of the single candidate event 
is not one of the most likely background channels, it is justified to give a cross-section for this process. 
Based on the single event, the production cross-section relative to CC interactions is 1.95^ 3 44(stat) ± 
0.29(syst) x 10~ 4 with a systematic error of 15% from the flux normalization. The rate of CC associated 
charm production with the threshold cut is obtained to be (4.50I3 44 (stat) ^ 0.68(syst)) x 10 -4 . The 
cross-section predicted by the QCD inspired parton model [ 5 ] has a strong energy dependence. Although 
the cross-section of this process at the average energy of the CHORUS beam is low, the measured 
production rate is in agreement with the prediction of the QCD inspired parton model. The relative rate 
of CC associated charm production is calculated to be ~ 2 x 10~ 4 within the framework of this model. 
The a posteriori probability that the background for the topology with two neutral decays gives one event 
is 0.016. Taking into account that this was not the only topology searched for, it is difficult to convert this 
into a uniquely defined confidence level. However, given the special topology, it is very likely that this 
event constitutes an observation of associated charm production in CC interactions. 
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